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Abstract 

We construct a particular class of quantum states for a massless, minimally coupled free 
scalar field which are of the form of a superposition of the vacuum and multi-mode two- 
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can produce an arbitrarily large amount of negative energy in a given region of space at a 
fixed time. This class of states thus provides an explicit counterexample to the existence 
of a spatially averaged quantum inequality in four-dimensional spacetime. 
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I. INTRODUCTION 



It is well known that quantum fields can produce local renormalized negative en- 
ergy densities. This is despite the fact that the classical expression for the energy 
density appears to be positive definite. The negative energy density is possible be- 
cause renormalization involves an infinite subtraction, and is needed to make the 
stress tensor operator well-defined. If there were no restrictions on the sorts of neg- 
ative energy densities attainable, one could have a number of bizarre possibilities, 
including traversable wormholes [|I|, 0, faster-than-light travel 0, J|, [f| time travel 
(I], |2], [5], and violations of the second law of thermodynamics 0. Such phenomena 
are, at best, rare. 

Fortunately, there are known to be severe constraints on the magnitude and dura- 
tion of negative energy fluxes or densities, at least in the case of linear fields. These 
are the quantum inequalities, which take the typical form || ||, [10|, ll|, [12, [13], [14]] 



/oo 
(:T 00 (x,t) :)F{t)dt>- 
-oo 



d 



Here (: T o(x, t) :) is the renormalized energy density in an arbitrary quantum state, 
F(t) is a normalized sampling function with a characteristic width to, c is a dimension- 
less constant, usually somewhat less than unity, and d is the dimension of spacetime. 
Most of the work on quantum inequalities is for the cases of d — 2 and d = 4. 
More precisely, Eq. (]!]) is a worldline quantum inequality for an inertial observer in 
Minkowski spacetime, where t is the observer's proper time. (Recently, worldline 
inequalities have been established in much greater generality in globally hyperbolic 
spacetimes The physical interpretation of this constraint is that there is an 



inverse relation between the magnitude of the negative energy and its duration. For 
example, in four dimensions, the observer will not see the negative energy last for a 
time longer than about \pmax\~ 1 , where p m ax is the maximally negative energy density 
seen by this observer. This type of constraint greatly limits the macroscopic effects 
of negative energy. In particular, macroscopic wormholes or "warp drive" spacetimes 
are severely constrained |16|, [H], pL8|| . Furthermore, one can use the worldline inequal- 
ities to place some constraints on the spatial distribution of negative energy ||19|| . All 



allowed spatial distributions must be such that Eq. (|l]) is satisfied for every inertial 
observer in Minkowski spacetime. 

There is a different type of constraint on the energy density, which is that the 
total energy in any quantum state for a quantum field in boundary-free Minkowski 
spacetime must be non-negative. Let the Hamiltonian operator be 

H = J d d - l x :T 00 (x,t) :, (2) 

where the integral is taken over all space at a fixed time. The expectation value of 
H in an arbitrary quantum state is non-negative: 

(H)>0, (3) 
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and (H) = only in the vacuum state. Thus any non-vacuum state must have net 
positive energy, even though there can be local regions with negative energy density. 

This raises the question of whether there are spatial or spacetime quantum in- 
equalities analogous to Eq. ([!]), but involving sampling over space or space and time. 
For spacetime averages, there is an abstract existence proof for a very general class 
of such inequalities in globally hyperbolic spacetimes. These include, for example, 



inequalities bounding violations of the dominant energy condition p0[ . However, at 
present quantitative estimates for these are not known. 

One can get some spacetime quantum inequalities by simply averaging the world- 
line inequalities over families of worldlines with spatial weighting factors, but in at 
least some cases these are known not to be sharp. For example, stronger constraints 
were found in two-dimensional Minkowksi space [JTTJ, ^TJ. We may express these as 



J dx dt F s {x) F T (t) (: T 00 {x, t) :) > -B(x , t ) . (4) 

Here Fg(x) is a spatial sampling function with width xq, Fr(t) is a temporal sampling 
function with width t , and the bound B(x ,t ) is finite as long as at least one of 
Xq and to is nonzero. The sampling functions are assumed to be each normalized to 
unity: 

J dxFs(x) = J dtF T (t) =1. (5) 

If the sampling functions are both Lorentzian functions, then B oc (xq + t )~ 2 ; if 
they are both Gaussians, then B oc (xq + tfy~ l . In either case, one obtains the 
corresponding worldline inequality when Xq = 0, and a nontrivial bound on the 
spatially sampled energy density when to = 0. One consequence of this is that in two 
dimensions the spacetime averaged results are stronger than the temporally averaged 
ones, since as to — > (for xo > 0) one gets a finite bound. 

It is natural to ask if there is a generalization of Eq. to four- dimensional 
spacetime, especially one which gives a nontrivial bound when the sampling is over 
space only. One of us |22| has given an argument that this is not the case, and 



that the spatially sampled energy density can be unbounded below. Actually, the 
argument was written out there for a more generic situation, that of evolution from 
a Cauchy surface in a general curved spacetime. While explicit constructions were 
given in this case, they were by pseudo differential operator techniques, and they were 
not translated into more conventional quantum field-theoretic terms. Also, the details 
of the more special case of a surface of constant time in Minkowski space were not 
written out. So no explicit version of the argument in conventional special- relativistic 
quantum field-theoretic terms has yet been given. The purpose of the present paper 
is to provide such an example and to draw as many physical insights as possible from 
it. In the following section, we will provide the explicit construction of a class of 
quantum states for the massless, minimally coupled scalar field, with negative energy 
density, as well as give the two-point function and the energy density in this class of 
states. We then show that although this state satisfies Eq. (JTj) , the spatially sampled 
energy density can be arbitrarily negative. Further implications of this example are 
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discussed in Sect. |TTT| . We work in units where h = c = 1. (A related construction has 
subsequently been used, in Ref . p3[ , to prove that there are no quantum inequalities 
along null geodesies in four-dimensional Minkowski spacetime.) 



II. THE ENERGY DENSITY IN A CLASS OF QUANTUM STATES 



A. Characterization of the State 

The particular quantum states which will be used in this paper are of the form of 
a superposition of the vacuum and a multi-mode two-particle state: 



\ij) = N 



|0)+ J2 Ck.kjkx,^) 

ki,k2 



(6) 



where we assume that the Ck 1 ,k 2 ' s are symmetric. Here |ki, k 2 ) is a two-particle state 
containing a pair of particles with momenta ki and k2, and the normalization factor 
is 



N 



(7) 



1 J2 l c ki,k 2 | 

ki,k2 

where £ = 2 if ki 7^ k2 and £ = 1 if ki = k2. Note that each of the two-particle 
states in Eq. (@) appears twice, since |ki,k 2 ) = |k 2 ,ki); this leads to the factor of 2 
in Eq. (|7|) when ki 7^ k 2 . We can quantize the scalar field in a box of volume V, so 
that the momenta are discrete. Later, we will let V — > 00, so that 



E 

k 



v 



(2vr) : 



d 3 k 



(8) 



In this limit, we replace the coefficients Ck l5 k 2 by functions of continuous momenta 
defined by 

V 

6(ki,k 2 ) = 77— r T7C kljk2 . (9) 



(27T) 



The normalization factor then becomes 



N 



1 + 2 / d^d 3 /^ |6(ki,k 2 ) 



(10) 



The quantum state is also characterized by a cutoff A, so that |ki|, |k 2 | < A. Later, 
we will consider the limit in which the cutoff is removed, A —>■ 00. We will require 
that the functional form of &(k 1; k 2 ) be such that the integral in Eq. (||) remain finite 
in this limit. 

The particular choice of 6(k 1; k 2 ) which we adopt is the following: 



6(k 1 ,k 2 )= x (k 1 + k 2 )(|k 1 ||k 2 |)^2 
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where 

( v J Xo, if |p| <Po f . 

X(P) = S (12) 
I , otherwise 

where Xo an d Po are arbitrary constants. (The function x(p) is taken to have a sharp 
cutoff merely for convenience; it could be made smooth without affecting our essential 
argument.) Thus our state is really a four-parameter family of states defined by the 
parameters A, u, xo, and Po- We will be especially interested in the limit in which 
the magnitudes of the momenta of both particles become large with p fixed. In this 
limit, the pairs of particles have almost exactly opposite momenta. 
The integral in Eq. (j5|) can now be written as 

Xlfd'kf d 3 p(|k||k-p|) 2 -\ (13) 

where k = kx and p = \a\ + k 2 . In the limit that A — ► oo, the k- integration will 
converge at the upper limit provided that 

v<-\. (14) 
This is the condition that our quantum state be normalizable in the A — ► oo limit. 



B. The Two-Point Function 



Now we wish to calculate the form of the two-point function for the massless scalar 
field in the states we have selected. First expand the field operator in terms of mode 
functions and creation and annihilation operators as 



(15) 



where the mode functions are 



i(k-x— u> t) 



with <jj — |k|. 

We are interested in the renormalized two-point function in the state \il 

by 

(: <p(x) V (x') :) = (iPMx) V (x'm - (0\<p(x) V (x')\0) . 
It may be expressed as 

(: V (x) V (x') :) = E[/kW/k'(^) (4«k') +Ux)f*,(x') (4,a k ) 



k.k' 



+ /k(^)/ k '(x') (a k a k ,> + fZ(x)fZ,(x') (44,) 



(16) 

defined 
(17) 
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where the expectation values are in the state \ip). 

If we evaluate this expression explicitly using the form of the quantum state in 
Eq. (|6|) and of the mode functions in Eq. (|16D, and then take the infinite volume limit, 
using Eq. (Rj), the result is 



2N 2 

-7T ruj / U K U K — -p 

(2vr)3 J yfav 
2e i(k'.x'-k.x) e i(a ) t-^'t') J d3jfe 1 6*(k 1 ,k)6(k 1 ,k') 



(: <p(x) V (x') :)= Rej d 3 k d 3 k' ^_ 



x 

+ e i(k.x+k'.x') e -i( W *Wt') 6 ( kjk ')j . ( 19 ) 

The first term on the righthand side, which is quadratic in b, comes from the expec- 
tation value in the two-particle component of alone, whereas the second term, 
which is linear in b, comes from a cross term involving both the vacuum and the 
two-particle component. It is the latter part which will be of greater interest to us. 
So long as A < oo, the integrals in Eq. fll9|) exist and define smooth functions of x and 
x', so the two-point function has the Hadamard form. More precisely, the behavior 
as x — > x' has the Hadamard form. The falloff as the separation between x and x' 
increases will be slow, on account of the sharp cutoff in x(j?)- Had we used a modified 
X, falling off smoothly, we could have arranged for a more rapid falloff as the spatial 
separation of x and x' becomes large. 



C. The Energy Density 

The energy density in a given quantum state may be obtained from the renor- 
malized two-point function. For the massless, minimally coupled scalar field, we can 
write 

p = \ lim [{d t d t , + V • V) (: <p(x) <p(x') :)] . (20) 

We may now use Eq. © to write the energy density as 
N 2 

~ (>o 3 

_ e i(k+k')-x e -i(a,+a/)i 6 ( kj k /)l ^ ( 21 ) 



P 



Re J rf 3 A;rf 3 A; , v / ^7(l + k-k / ) ^e^'^-V^')* J d 3 k±b* (k 1; k)&(k 1; k') 



where k and k are unit vectors in the directions of k and k, respectively. Let pi 
denote the term quadratic in b, and p2 denote the term linear in b. 

Let us first examine the behavior of p2- If we use the explicit form for b given in 
Eq. ([11]), and let p = k + k, we find 

Re[d 3 k I rf 3 pvW(l + k-k)e ip - x e- i(w+|p - k|)< (w|p-k|)^^. 

J J\p\<Pa 

(22) 



XoN 2 

P2 



(27T) 
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In the limit that uj becomes large compared to po, we can expand the product k • k 
as follows: 



k-k' = , k(p ~ k) =(-i + ±v 



oj 



k • p p 2 
1-2 — - + ^ 

UJ UJ Z 



(p-k)-(p-k) 

p 2 -(k-p) 2 , . 

To leading order when 3> po ; we can write uj' ~ u;. If the integral for p2 is dominated 
by values of uj large compared to po, we have 



XoN 

P2 



2 



, Re / rf 3 A; / d 3 p uj^-V [p 2 - (k • p) 2 ] e ip ' x e" 2 ^' . (24) 
2{2n)- i J ^IpI^po 

We next perform the angular integrations. If we integrate over the directions of k 
with p fixed, we find 

J d 3 kp 2 = 4np 2 J dujuj 2 (25) 

and ^ 

J d 3 k(k-p) 2 = 2np 2 J duuj 2 J dec 2 = yp 2 J ' dujuj 2 , (26) 

where c = cos# and 9 is the angle between k and p. Next we integrate over the 
directions of p with x fixed to find 

f 2 (po,r) = [ d 3 pp 2 e ipx = 2n [ P dpp A [ dc' e iprc ' = — [ P dpp 3 sin(pr) 
j|p|<po Jo J—i r Jo 

= ^[3(plr 2 ~2)sm(p r) -p r(plr 2 -6)cos(p r)}, (27) 

where r = |x|, d = cos a, and a is the angle between x and p. The function f 2 (po, r ) 
has the following behavior: 

-(4vr/r 2 ) p 3 cos p r, r > po -1 , ooX 
fc(Po,r) ~ { ■ (28) 

(4tt/5) pi , r < po 1 



Now we can write Eq. (|24|) as 



YnN 2 f A 

P2«-^-/ 2 (po,r) / dujuj 2 " cos(2ujt), (29) 

07T J q 2 po 



where q 2 is a constant chosen so that the approximations made in finding the integrand 
(uj ^> po) are valid throughout the range of integration. Note that if t — and u > — |, 
the integral will diverge in the A — > oo limit. 
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Let us now set t = and v = —1/2. Then the cu-integral becomes 

/ — = In , (30) 



lq 2P0 U \q2P0 

so that 



P2«-^/ 2 (po,r)lnf— ) . (31) 

We now examine the contribution of p\. From Eq. fl2"l|) , we have 
AT2 



Pi 



(27T) 



Re y Aft'vW(l+k-k') J 2e i(k'-k).x e i(^')t y d 3 A; 1 &*(k 1 ,k)6(k 1 ,k / ) 



(32) 

If we use Eqs. ( ITT] ) and (|T2]), letting p = k x + k and p' = k x + k', and interchanging 
the order of integration, we can write 

Pi = Re f d% {k^- 1 [ d 3 p [ d 3 p> (1 + k ■ k') 

(ZTT) 6 J J\P\<P0 Ap'\<P0 

x (Ip-kx^lp'-^D^e-^-P'^e^P-^l)*-^'^ 1 ^. (33) 

If the kx-integral is dominated by values of | ki | ^> Po, then we can use |p — ki | v w 
|p' — kil^ « |kiT- Therefore, at t = we have 

Pi- 2 ^° 2 ^ 2 Re I d*k x \k x \ Av - 1 j d 3 p [ d 3 p' (l + k-k')e- i(p - p ' )x . (34) 
(2?r) 3 7 J|p|<po Vl<po 

Note that in the high k\ limit, 

k . k ' = (p - ki) • (p' - ki) ^ 

'(p_ kl ).(p_ kl ) /(p/_ kl ).(p/_ kl ) 



As before, we now set v = —1/2. A similar calculation to that given previously 
for p 2 yields 

^ 2 Xo 2 iV 2 / A \ 

7T \qm) 

where q\ is a constant defined analogously to q 2 , and 

/i(Po, r ) = — — [sin(p r) - Po?" cos(p r)] 2 , (37) 

and 

(47rp ) 2 / r4cos2 (Por) , r > PtT 1 , 00 . 
fi[Po,r) ~ ^ . (38) 

47r/3 2 p6 r<p ^ 
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Both |/0 2 1 an d pi diverge logarithmically, but we can arrange for \p2\l pi > 1, provided 
that 




where 



Y 



In (A/ q 2 p ) 



as A • x . (40) 

_ln(A/g!p )J 

This is possible because p\ and p 2 involve different powers of Xo, stemming from the 
fact that while p 1 is quadratic in b, p 2 is only linear in b. (We note in passing that 
the choice of v — —1/2 gives xo dimensions of length.) 

For any value of po, the energy density is approximately constant in space over a 
region of size, r <C po^ 1 . In this region, (/2//1) ~ 9/(20 irp ). Choose 

Xo < 1 (£W « . (41) 

where we have used the fact that Y — > 1, as A — > 00. This allows P2 to dominate 
Pi, so that p = pi + p 2 can be made arbitrarily negative over an arbitrarily large 
region. In summary, we have a normalizable state in which the local energy density 
may be made arbitrarily negative at time t = in a finite region of space. It should 
be pointed out that in the strict A — > 00 limit, our state will no longer have the 
Hadamard form. However, the key point here is that we can consider a sequence of 
states, each with a finite value of A. Therefore each state in this sequence will have 
the Hadamard form. By progressively increasing the values of A as we vary over 
the states in the sequence, we can make the energy density in our spatially sampled 
region as negative as we like. In the next section, we will discuss some of the insights 
which may be drawn from this example. 



III. IMPLICATIONS AND DISCUSSION 

We are now in a position to draw a number of conclusions from our example. The 
first conclusion is the proof that there are no lower bounds on the spatially sam- 
pled energy density for the massless minimally coupled scalar field in four spacetime 
dimensions. To be more precise, consider the spatial average of the energy density 

S(F) = |rf 3 xF(x)(:T 00 (x,t):) (42) 

at time t, where F(x) is a spatial sampling function. Let r be the characteristic 
width of this function. More generally, tq can be the maximum width of F(x) in any 
direction. If there were to be spatial quantum inequalities here, S(F) would have to 
have a lower bound when we vary over all quantum states for a given F(x). However, 
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the energy density in the above example varies in space only over a scale of the order 
of 

Ac = - , (43) 
Po 

which can be arbitrary. Furthermore, we can make the energy density in a region 
whose size is less than Ao arbitrarily negative. For a given sampling function F(x), 
we can choose Ao > r , so that the energy density is approximately constant over the 
region being sampled and the spatial sampling has no effect, yet the energy density 
in that region is arbitrarily negative. Thus S(F) is unbounded below. However, this 
state satisfies the temporal quantum inequality, Eq. (|l|). One can argue from energy 
conservation that physically what is happening is that there must be large fluxes of 
positive energy entering the region, which damp out the negative energy sufficiently 
for the temporal inequality to hold. Recall from Eqs. (^) and ([33|) that p 2 oscillates 
much more rapidly in time than does p\. Thus if p 2 < at one time, it has changed 
sign a short time later. This causes the time integral of p 2 to tend to average to zero, 
and allows the worldline quantum inequalities to be satisfied. 

We can gain more physical insight into the behavior of this class of states by 
recalling that in the high frequency limit, the pairs of particles in the two-particle 
component of the state have nearly but not exactly opposite momenta, as shown in 
Fig. |T[ This allows the scale for the spatial variation of p to be large compared to 
the temporal scale, which is set by the characteristic frequency of the excited modes. 
This can be seen in Eq. fl^f), where the spatial factor e* px oscillates much more 
slowly than the temporal factor Q- 2luJt . One can summarize the difference between 
the spatial and temporal behaviors by noting that it is possible for the three-momenta 
of a pair of particles to cancel, but not their energies. We can also see that it is a 
very special property of two-dimensional spacetime which allows for the existence 
of spatial bounds. Here the momenta of a pair of particles must be either exactly 
parallel or exactly anti-parallel. However, if they are anti-parallel, there is a factor 
analogous to the factor of 1 + k • k in Eq. (|21~D which suppresses their contribution. 
Although in this paper we have looked only at four- dimensional spacetime, it seems 
likely that one can construct similar counterexamples to spatial quantum inequalities 
in all spacetime dimensions greater than two. 

We can also understand why the Hamiltonian, obtained by integrating over all 
space, is bounded below, yet the spatial average over any sampling function with 
finite width is not. The integration over all space removes the contribution of the p 2 
term, leaving only the non-negative contribution of the p\ term. Thus the net energy 
in any state must be non-negative. However, it is possible for the spatial distribution 
at a fixed time to be such that the compensating positive energy is arbitrarily far 
from the negative energy, and hence violates a spatial quantum inequality. 

The current results also shed light on an earlier proposed explanation of the 
"Garfinkle box", given in Ref. [I9|. This refers to an unpublished result of Garfin- 
kle pi, who showed that the total energy of a scalar field, E, contained within an 



imaginary box in Minkowski spacetime, at fixed time, is unbounded below. The box 
is "imaginary" in the sense that there are no physical boundaries at the walls of the 
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p = k + k' 



k 



k' 



FIG. 1: The pairs of particles in the two-particle component of the state have nearly but 
not exactly opposite momenta, in the high frequency limit. Although their momenta nearly 
cancel, their energies do not. 

box. The Garfinkle result is that there exist quantum states for which E is arbitrar- 
ily negative. The proposed explanation given in Ref. [19j for the unboundedness of 
the total energy E in this box involved two factors: (1) The energy is measured at 
a precise instant in time, and (2) the walls of the box are sharply denned. It was 
argued that this allows an arbitrary amount of negative energy to have entered the 
box by time t, while at this time excluding an even larger amount of positive energy 
which may be just outside of the box at time t. The results of the present paper 
show that, although this is an explanation of the Garfinkle result, it is not the only 
explanation. In particular, the unboundedness of the energy in our case does not 
require the spatial region over which we average to be sharply defined. 

In this connection, we emphasize that our results show that the unboundedly 
negative energies cannot be ascribed to "edge effects" associated with the sampling 
function. This is because in our examples, the energy density is uniformly negative 
throughout the sampling region. 

While spatially averaged quantum inequalities do not exist, we noted earlier that 
spacetime averaged ones do. This means that the lower bounds associated with the 
spacetime averages must diverge as the temporal averaging scale goes to zero. It is 
natural to conjecture that this is related to the divergence in the spatially averaged 
inequalities, as follows. Suppose the most severe divergence attainable in the spatial 
average is 

J dx Fs(x) (: T 00 (x, t) :) diverges as —/(A) as A — > oo , (44) 

for some function / of the ultraviolet cutoff A. Similarly, suppose that for a temporal 
averaging function Fx(t) with characteristic width t , we have 

J dx dt Fs(x) Fxif) (: T o(x, t) :) diverges as — g(t ) as to — ► , (45) 

for some function g(to). It is natural to ask if /(A) ~ g(l/A). A relation like this 
would probably give the most direct physical significance to the divergences in the 
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spatial averages. (Such a relation would not hold in two dimensions, but might in 
four.) 

The logarithmic divergences (of the energy density with the cutoff A) are rather 
slow. We do not know whether it is possible to construct examples in Minkowski 
space with more severe divergences. (In Ref. [2~2"|, more severe divergences were 



found at Cauchy surfaces with non-zero second fundamental forms. In generic curved 
spacetimes, this would be the case for all Cauchy surfaces.) 

In summary, we have shown that there are no purely spatially averaged quantum 
inequalities over bounded regions in four- dimensional Minkowski spacetime, even 
though the integral over all space is bounded. However, spacetime averaged quantum 
inequalities do exist. As noted earlier, in at least the two-dimensional case, spacetime 
averaging seems to lead to tighter bounds than does temporal averaging alone. Recent 



numerical results of Dawson and Fewster 251 indicate that similar behavior occurs 



in four dimensions. However, at present the extent to which spacetime averaging 
improves temporal averaging is unclear. Thus there appear to be two approaches for 
better understanding the allowed spatial distributions of negative energy. One is the 
search for spacetime averaged quantum inequalities, and the other is the continuation 
of the program begun in Ref. fl9f , which seeks to extract as much information as 



possible from the worldline inequalities. Both of these approaches are currently under 
investigation. 
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